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Abstract 

We study the average number A n per site of the number of differ- 
ent configurations of a branched polymer of n bonds on the Given- 
Mandelbrot family of fractals using exact real-space renormalization. 
Different members of the family are characterized by an integer pa- 
rameter b, 2 < b < oo. The fractal dimension varies from log 2 3 to 2 
as b is varied from 2 to oo. We find that for all b > 3, A n varies as 
\ n exp(bri^) : where A and b are some constants, and < ip < 1. We 
determine the exponent ip, and the size exponent v (average diameter 
of polymer varies as n u ), exactly for all b, 3 < b < oo. This generalizes 
the earlier results of Knezevic and Vannimenus for 6 = 3 [Phys. Rev 
B 35 (1987) 4988]. 

1 Introduction 

The study of statistical physics models on deterministic fractals has a long 
history [1, 2, 3, 4]. Linear and branched polymers on fractals with finite 
ramification number provide very simple and pedagogical examples of renor- 
malization group techniques at work: these system show a non-trivial critical 
point, and the values of the critical exponents can be determined by lineariz- 
ing the exact real-space renormalization transformation. The renormaliza- 
tion equations are coupled polynomial recursion equations in a finite number 
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of variables, and are easy to study. By studying different geometrical fractals, 
one can investigate how the critical exponents change with the geometrical 
properties of the underlying space. 

One particular family of fractals which has been used often for such stud- 
ies is the Given-Mandelbrot family of fractals [5]. Different members of the 
family are characterized by an integer 6, with 2 < 6 < oo. As we increase 
b from 2 to oo, the fractal dimension increases from log^i to 2. The criti- 
cal properties of linear polymers on the b = 2 fractal were first studied in 
[2], and these results were extended to b < 8 by Elezovic et al [6] using the 
exact renormalization equations. Surprizingly, it was found that while the 
exponent appeared to converge to the two-dimensional value 3/4, as 6 was 
increased from 2 to 8, the difference in the susceptibility exponent 7& from 
the known exact value 43/32 in two dimensions was found to increase with 
increasing b. This was explained in [7], where the asymptotic behavior of 
critical exponents for large b was determined theoretically using finite-size 
scaling arguments, and it was shown that 7& should tend to a different value 
133/32 for large b . Numerical Monte Carlo renormalization group techniques 
have been used to estimate the critical exponents for significantly larger val- 
ues of b up to 80 [8, 9]. Knezevic and Vannimenus (KV) used the real-space 
renormalization technique to study the properties of branched polymers on 
the 6 = 2 fractal, and also studied the transition from the extended phase to 
collapsed phase [10]. This was later extended to other fractals, including the 
6 = 3 fractal [11]. Dense branched polymers for the 6 = 2 fractal have been 
studied in the context of spanning trees and loop-erased random walks [12], 
and the abelian sandpile model [13]. However, a study of branched polymers 
on fractals for higher 6 has not been undertaken so far. Nor are the properties 
of the large-6 limit known. 

In this paper we study the number of different configurations of an n- 
bond branched polymer on the Given-Mandelbrot family of fractals using the 
exact real-space renormalization group techniques. On regular lattices, this 
number usually varies as \ n n~ e , where A is some lattice -dependent constant, 
and 9 is a critical exponent. General theoretical arguments that prove the 
exponential growth would allow stronger correction terms like exp(6n^), with 
ip < 1. Why the first correction term to the exponential growth is a simple 
power-law term is not fully understood. To see how general is the power-law 
correction form, one can study this question on different graphs, e.g. fractals. 
We find the power-law correction also on the 6 = 2 fractal. However, this 
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case is exceptional. For all 6 ^ 2 , while the number of configurations still 
increases exponentially with n, the leading correction term to the exponential 
growth is the stretched-exponential form: this number varies as A"e bn * , where 
A and 6 are some constants, and < ip < 1. We determine the singularity 
exponent ip, and the size exponent v (average diameter of polymer varies as 
n u ), exactly for all b, 3 < b < oo. This generalizes the earlier results of KV 
for 6 = 2 and 3. 

This paper is organized as follows: In section 2, we start by recapitulating 
the definition of the Given-Mandelbrot family of fractals, and introduce the 
generating function for the number of branched polymer configurations with 
n monomers. Since the fractal does not have translational invariance, we 
average over different positions of the polymer. The general technique of 
real-space renormalization applied to these problems is outlined in Section 3, 
using the 6 = 2 case as an illustrative example. The qualitative behavior of 
the renormalization equations for b > 3 is discussed in Section 4. It turns out 
that while the equations involve rather complicated high-degree polynomials, 
the critical exponents v and ip do not depend on most of the terms in these 
polynomials. We can ignore most of these terms, and still determine the 
exact values of these exponents, if we can identify the "dominant terms" in 
the recursion equations. This is done in Section 5. Finally, in Section 6, 
using our knowledge of the dominant terms, we determine the exponents v 
and ip for all b > 3, without having to write down the full set of recursion 
equations. 

2 Definitions 

For any given integer b, 2 < b < oo, the recursive construction of the Given- 
Mandelbrot family of fractals is shown in Fig. 1. We start with a graph 
with three vertices and three edges forming a triangle. This is called the 
first order triangle. To construct the graph of the (r + l)-th order triangle, 
we take graphs of 6(6 + l)/2 triangles of r-th order, and glue them together 
(i.e. identify corner vertices) as shown in the figure, to form an equilateral 
triangle with base which is b times longer. The case 6 = 2 corresponds to the 
well-known Sierpinski gasket [Fig. 2]. 

It is easy to see that the number of edges in the graph of the r-th order 
triangle is 
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Figure 1: The recursive construction of the Given- Mandelbrot fractal for 
6 = 4. (a) The graph of first order riangle.(b) the graph of a (r + 1) order 
triangle, formed by joining b{b + l)/2 r-th order triangles shown as shaded 
triangles here. 




Figure 2: The graph of a 5-th order triangle for b = 2. 
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3b r (b + l) r 2~ r , and all vertices have coordination number 4 or 6, except 
the corner vertices. The distance between the corner vertices of r-th order 
triangle is b r ~ l . Thus, the fractal dimension of the graph is D b = \og b [b(b + 
l)/2]. For b = 2,3,4.., these values are 1.5849, 1.6309, 1.6609... respectively. 
For large b, the fractal dimension tends to 2 as D b pa 2 — log b 2. The spectral 
dimension D b of the graph can also be calculated exactly for general b [14]. 
The values of D b for b = 2 to 10 are listed in [15]. For large b, D b tends to 2, 
and the leading correction to its limiting value is given by D b pa 2 — lo f lo f b 

[16]- 

The determination of the generating function for the branched polymers 
on these fractals follows the treatment of [2, 10, 11]. Let A n (N) be the 
number of distinct single connected cluster of n bonds on a graph with N 
total number of bonds in the graph, different translations of the cluster being 
counted as distinct. For large N, this number increases linearly with N. We 
then define 

A n = IAm N ^ 00 A n (N)/N; (1) 

and 

oo 

G(x) = ]T A n x n , (2) 

n=0 

where we assume the convention A — 1. 



3 Renormalization equations for the 6 = 2 
fractal 

We assign a weight x n to each configuration of the polymer with n occupied 
bonds, and define restricted partition functions A^ r \ B^'\ C^'\ D^ r \ and 
F^'\ as shown in Fig 3. Here A^ is the sum of weights of all connected con- 
figurations of the branched polymer inside a r-th order triangle, such that 
only one of the corner vertices is occupied by the polymer, and the other two 
corner vertices are unoccupied. is the sum of weights of all configura- 
tions with two mutually disconnected clusters, each cluster connected to a 
specified corner vertex. consists of sum of all configurations of polymer 
that connect two specified corner vertices of the triangle, with the third ver- 
tex remaining unoccupied. Similarly, an F^ correspond to sums over 
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Figure 3: Definition of the restricted partition functions 

A<r\ fiW, C< r \ £> (r) , £W anc i fW. 



all configurations with all three corner vertices occupied, and mutually dis- 
connected, and connected respectively. corresponds to all three corner 
vertices occupied, but only two specified vertices connected to each other by 
paths involving occupied bonds lying within the triangle. 
The values of these restricted functions for r = 1 are 

A W = D W = £(i) = i; £(i)= C^=x; F^ = 3x 2 + x 3 . (3) 

It is straightforward to write down the recursion equations for these re- 
stricted partition functions at level (r + 1) in terms of those at level r. One 
has to sum over all possible polymer configurations for different r-th order 
triangles, subject to the constraint that all occupied bond are connected to 
one of the corner vertices of the (r + l)-th order triangle. For example, for 
b = 2, KV obtained 
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A (r+i) = A][l + 2 B + 2B 2 } + 2B 2 C + F[B 2 + A 2 + 2BD] (4) 
B (r+i) = B 2 + B 3 + p^ ABC + 2AB] + F 2 [B + £>] (5) 

^(r+1) = AB i + 352^ + F [ 7C <2 + + + £j ^ 

£,(r+i) = A 2 + B[6C 2 +4AC+2A 2 ]+D[2B + W 2 ]+2F[2CD+AD + BC+BE] 

(7) 

E (r+i) = ^3 + UC 3 + 12j bc7) + 6ABD + 3B 2 £ + 3F[C 2 + D 2 + 4C£] (8) 

F (r+1) = 3F5 2 + g F 2 c + F 3 (g) 

where we have dropped the superscript (r) over all the variables in the right- 
hand side of all the recursion equations to simplify notation. The generating 
function G(x) is expressible in terms of these variables 

oo 

G(x) = ]T 3- r (A^ 2 + AW'bW + B^D^) (10) 

r=l 

It was shown by KV that there exist a critical value x c , such that for all 
x < x c , (A( r \ fiW, C (r \D^ r \ E^ r \ FW) tends to a fixed point (A*(x), 0, 0, A*\x), A*\x), 0), 
where the value of A*(x) increases monotonically from 1 to oo as x increases 
from to x c . For all x > x c , A<r\ B<r\ C<r\ D<r\ E<r\ all diverge to 
infinity. At x — x c , the values of A^ r \C^ r \D^ and E^ diverge to infinity 
for large r. If we change variables to A = AF, B = B,C = CF, D = DF 2 
and E = EF 3 , F = F, the non-trivial fixed point occurs at finite values of 
A, B, C, D, E, with F — 0. Linearized analysis of the renormalization equa- 
tions near the fixed point determines the singularity exponent 6 for the func- 
tion G(x) ~ Ixc — x] 6 ^ 1 , and the exponent v ( correlation length ~ \x c — x\~ u ). 
KV found v ~ 0.71655, and 6 ~ 0.5328. 

4 Renormalization equations for b > 3 

Interestingly, the qualitative behavior of the recursion equations is very dif- 
ferent for b > 2. For 6 = 3, the singularity of G(x) is not a power-law 
singularity, but an essential singularity [11]. The case b > 3 has not been 
studied so far. 

For a general value of b, the equations are still coupled polynomial recur- 
sion equations in six variables. We denote the six functions A^'\ B^'\ . . . F^ 
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by K\ , with i — 1 to 6, and represent the polynomial recursion equations 
schematically as 

Kt +1) = m^ r) }), for i = 1 to 6. (11) 

where fi are some (6-dependent ) polynomial functions of their six arguments. 
One can also write the generating function G(x) in terms of {K^} in a 
polynomial form similar to Eq.(lO). 

oo 

G(x) = E [b(b+ l)/2r g({K^}) (12) 

r=l 

It is rather tedious to write down the explicit forms of these polynomials 
for any 6 > 3. The number of terms in each increases very fast with 
6. There are approximately 100 terms in each of the recursion equations 
for b = 3 [17] , and the number would run into thousands for b = 4. The 
number of terms would increase as b 12 for large b, as that is the number of 
polynomials with six variables with maximum degree b(b + l)/2. Also, the 
coefficients of the terms become very large, increasing as exp(b 2 ). Even for 
b = 3, to generate the recursion equations, one has to use a computer. It 
seems clear that a brute-force approach to determine the recursion equations 
is not feasible, except for a few additional values of b. 

Interestingly, even though the recursion equations are rather complicated, 
we show below that for all b > 3, the generating function has an essential 
singularity of the type G(x) ~ exp(a\x c — x\~ a ). This corresponds to A n ~ 
x~ n exp(bn^), where b is some constant, and ip = a/(l + a). We determine 
the exact value of the 6-dependent exponents ip and v for all b. 

If we start with the initial conditions (Eq. (3)), and iterate the recursion 
equations, KV found that many qualitative features of the behavior of the 
recursion equations for the 6 = 3 case are same as for 6 = 2: for all x 
below a critical value x c , (A^ r \ B^\C^'\ D (r \ E (r \ fW) tends to a fixed 
point (A*(x), 0, 0, A* (x),A* (x),0), where A*(x) diverges to infinity as x 
tends to x c from below. For all x > x c , the values of all {Ki} tend to infinity 
for large r. For x = x c , both for b = 2 and 3, the values of A (r \ and 
tend to infinity, while and decrease to zero with iteration. The main 
difference is the limiting behavior of at x = x c . It tends to a non-zero 
limit for 6 = 2, but to zero for 6 = 3. Also, the variables {logK^} increase 
linearly with r for 6 = 2, but they increase exponentially with r for 6 = 3. 
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What makes this problem tractable is the fact that while the polynomial 
recursion equations are complicated, the asymptotic behavior of the variables 
depends only on a few terms in the recursion equations. To see this, consider 
first the simple case when each function /j has only a single term, and is of 
the form 

K? +1) = Ci {[[K^ (13) 
3=1 

where m is 6 x 6 matrix of non-negative integers. These recursion equations 
reduce to linear recursions on taking logarithms of both sides. We get 

\ogK\ r+1) = £ log iff + logci (14) 

3=1 

(r) 

These are easily solved. We get that for large r, the leading behavior of K\ 
is given by 

iog^ (r+1) = D(m% log + i^zrh lo s c il ( 15 ) 

3=1 

Let the eigenvalues of the matrix m be {X a }, a = 1 to 6, ordered so that 
X a > X a +i- Let the corresponding right-eigenvectors be Vi a , so that 

™.ijVi a = \ a V ia (16) 

Then, using the eigenvector decomposition of m, we see that for large r 

logK^ = 5iX[vn + 82X2^2 + higher order terms, (17) 

where Si and S2 are some coefficients, that can be expressed in terms of 
Cj's, and the left eigenvectors of m. As m has non- negative elements 
and Cj > l,vn are of the same sign, which can be chosen positive. Then all 
KiS increase or decrease with iteration for large r according as Si > or 
Si < 0. Thus Si — must correspond to x = x c , and for small deviations of 
x from x c , by continuity, we will have Si proportional to (x — x c ). 

If Si = 0, then the behavior of K^s for large r is governed by the second 
term in Eq.(17). In this case, all i^'s are not of same sign. If V12, v 42 and t> 52 
are positive, and the rest negative, we would get Ki,K± and K 5 to diverge, 
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and K 2 , K 3 and K 6 tend to zero for large r, as expected from the preceding 
discussion. 

Conversely, consider the full recursion equation(ll). For any two infinite 
sequences and whose logarithm diverges to infinity as r increases, we 
define the notation 

lo£f (7,( r ) 

a b, iff lim f^-T = !■ (18) 

We make the ansatz that at x = x c , for large r we have 

K\ r) = exp(vA r ). (19) 

Then, as the number of terms in the recursion equations is finite, for each 
Ki, there must be at least one term in the right hand side of its recursion 
equation for which the asymptotic rate of growth of the logarithm is exactly 
as the same as that of the left hand side. We can define a matrix m such 
that this dominant term is of the form of Eq.(13). Then we must have 

6 

\v i = ^2m ij v j , (20) 
i=i 

while all other terms in the equation for Ki are either neglible for large r, 
or make a contribution of comparable amount. These dominant terms define 
the matrix m. The vector {vi}, then , is a right eigenvector of m, with 
eigenvalue A. 

If there are more than one terms that make contributions of the same 
order, say in the recursion equation for Ki, then for any two such terms, 
there are two different row- vectors rrii and m- that satisfy eq.(20). Then we 
can subtract these to get a relation 

6 

EKr m 'yh = ' ( 21 ) 

There is one such equation for every such pair. However, all of them 
are not independent. We first reduce these to the minimal set of linearly 
independent equations, then, each one of such equations can be used to 
eliminate one of the i>'s, and then write Eq.(20) as an equation in fewer 
variables, with a new lower-dimensional matrix m. This also changes the 
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coefficients q and the eigenvectors, but does not change the value of Ai and 
A 2 . 

For example, for the 6 = 3 fractal, KV found that the dominant terms in 
the recursion equations are 



A' = 2A 2 B 2 


(22) 


B' = A 2 B 2 F + 3AB 4 


(23) 


C = 3A 2 B 4 + A 3 B 2 F 


(24) 


D' = AA 3 B 2 


(25) 


E' = 6A 4 B 2 


(26) 


F' = 2B 6 + QAB 4 F 


(27) 



We note that the right hand side involves only A, B and F, and the 
variable C,D and E can be determined in terms of these. In Eq.(23), the 
two terms on the right hand side are of same order iff AF = B 2 . This 
corresponds to the condition 

ui - 2v 2 + v 6 = (28) 

Note that we get the same condition if we had used Eqs.(24), or (27) instead 
ofEq.(23). Let 

Air) rp(r) 

lim - -. — = f* (29) 

Then we can write F = f*B 2 /A for large r, and eliminate F and work 
with a simpler set of recursions 

A' = 2A 2 B 2 , B' = (3 + f*)AB\ (30) 

This corresponds to a 2 x 2 matrix 

m=(H) < 31 > 

With this change of variables, the other equations also reduce to single term 
equations. Also, the matrix m, and hence the eigenvalues Ai and A2 do not 
depend on the coefficient (3 + /*). 
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Similarly, in other cases, the asymptotic behavior of the variables K^s 
with the full complicated polynomial recursion relations may be reduced to 
that for a simpler system of equations where only one term ( to be called the 
dominant term) is kept from each polynomial fa in Eq.(ll). 

The dominant terms are not unique. Only a few terms in each of the poly- 
nomials fi({Kj}) are dominant, and any of them can be used for determining 
the matrix m, and hence the eigenvalues Ai and A2. 

We note that in the neighborhood of different fixed points, different 
terms are dominant. For example, for branched polymers with attractive 
self-interaction, the fixed point corresponding to the dense phase is {K{\ = 
{0,0,oo,0,oo,oo}. Clearly, near this fixed point, the dominant terms are 
different. 



5 Identifying the dominant terms 

The problem of determining the critical exponents for our problem thus is 
reduced to that of identifying what are the dominant terms in the recursion 
equations. Each term in the recursion equation corresponds to a class of 
configurations of the polymers. However, all possible combinations of powers 
are not allowed in a given equation, as there are connectivity constraints on 

6(6+1) 

the allowed configurations. For example, a term like E 2 corresponding 
to all r-th order triangles with configurations of type E has monomers not 
connected to the corners, and is not allowed. 

Also, we note that while we may allow polymers with loops, the ratios like 
p{r) jQ{r) an j g(r) 1 £)(r) ^ en( \ ^ Q zero f or l ar g e r . Thus, any term correspond- 
ing to a polymer with loops is dominated by one without loop, obtained by 
removing one of the bonds in the loop ( which changes a type F triangle 
into type C, or type B into type D). Thus, loops are irrelevant, and the 
dominant terms correspond to polymer configurations without loops. 

To understand the characteristics of dominant configurations better, it is 
instructive to look at configurations that differ from each other locally. 

Consider configurations C\ and C2 shown in Fig. 4. We assume that these 
are identical to each other, except in the three r-th order triangles shown. 
In Ci, the vertex common to the three triangles shown is not connected to 
the polymer, but in C2 it is. Then if C\ is an allowed configuration, which is 
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connected to the outside in a specified way, so is Ci- If the weight of the rest 
of polymer outside the three triangles is W, the weight of C\ is WA, while 
that of C 2 is WBA 2 . The ratio of weights is AB, which tends to infinity for 
large r ( we shall show this later). This implies that a configuration of type 
C2 will dominate over the corresponding C±. 

Again, assume that C3 and C4 in Fig. 4 are the same configuration 
outside the three triangles shown. The ratio of weights of C 4 and C 3 is 
A 2 F/B. This also tends to infinity for large r (also proved later), and hence 
given a configuration of type C\ or C3, with two 'empty' triangles near the 
polymer, we can attach a sidebranch to the polymer to create a configuration 
that dominates over it. 

This technique of creating a more dominant configuration by attaching 
sidebranches works only if the two triangles into which the polymer is ex- 
tended are initially totally empty of polymer. This is shown in Fig. 5, where 
again we obtain the configuration Cq by attaching more monomers to the 
polymer configuration C5, with configuration of polymer outside the trian- 
gles remaining unchanged. In this case, the ratio of weights of Cq and C5 is 
DB/A, which we shall show tends to zero or large r. Hence configurations of 
type C5 dominate over C§ . Similarly, we can see that ratio of weights of Cg 
and C7 is C, which tends to zero for large r. Thus, a configuration like C7 
is dominant over the corresponding configuration of type Cg. We see that C 
or D type vertices are not favored in creating a dominant configuration. 

We can start with any allowed configuration of the polymer, and use 
this local modifications to generate configurations which are more dominant 
(Fig. 6). If we start with an initial configuration with a small segment of 
polymer, with many empty triangles, we find growth at tips or sidebranches 
is favorable. We continue this until no such growth sites can be found, and 
any further growth of polymer reduces its weight. This is then the maximal 
weight configuration. Fig. 7 shows two such configurations for the 6 = 7 
fractal corresponding to configurations of type A. We see that in a maximal 
weight configuration in a (r + l)-th order triangle, the polymer goes through 
as many as possible of the corner vertices of the r-th order triangles that 
are inside the (r + l)-th order triangle ( not at the boundary of it), and do 
not contain any type C, D or E vertices. This is a generalization of the 
observation of KV [11] that the dominant terms in the recursion equations 
had the central node of the triangle occupied. 

To generate the different terms that are dominant for the different terms 
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Figure 5: Some examples of configurations where extending the polymer 
configuration is not favorable. Here C5 has higher weight than Cq, and C7 
has higher weight than Cg. 




Figure 6: Generating the dominant configuration by extending the cluster. 
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Figure 7: Two of the dominant configurations that contribute to A^ r+1 ^ for 
the b = 7 fractal. 



in different Ki, we can take the dominant configuration for A^ r+1 \ and mod- 
ify it appropriately. For example, consider the dominant configurations for 
B^' +l \ Now, in the dominant configuration for A^ r+1 \ the polymer is as 
extended as possible, and hence will reach very close to the other corners of 
the (r+ l)-th order triangle. Then we have to make only a local modification 
in the triangles near this corner to connect it to the corner. This is shown in 
Fig. 8. Let the weight of the dominant configurations of A^ 1 ^ be X, and 
the multiplicative factor needed to convert its weight into that of I?( r+1 ) be 

y. 

Then, we have, by definition, 

A (r+1) s X . B (r+1) xy (32) 

To get the dominant configuration of F^ r+1 \ we have to connect the 
dominant polymer configuration of A^ r+1 ^ to both of the other corners. As 
the polymer in the dominant configuration reaches close to both of them, 
these local modifications can be done independently. Hence we get 

p(r+i) s X y 2 (33) 

To get the dominant configurations for D^ r+l \ we have to add a polymer 
segment disconnected to the first to the second corner. This can also be 
done by a local modification of the configuration near the corner ( see Fig. 
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8). We define Z as the factor by which we have to multiply the weight of 
the configuration of A^' +lS) to affect this change. Thus, we have 

D {r+i) XZ (34) 

Dominant confugurations for i^( r+1 ) involve two such local changes, which 
can be done independently. Similarly, for C (r+X) also, we have to make local 
modifications at the two corners, connecting one corner to the existing A^' +1 
cluster, and adding a disconnected cluster at the second corner. This gives 
us 

E (r+1) ^ XZ 2. C {r+1) ^ xyz ( 35 ) 

Thus, as far as the dominant terms are concerned, we can express the six 
functions f\ in terms of three functions X, y and Z. Equivalently, we can 
eliminante three of the variables (say C, E and F) using the relations 

D 2 = AE, B 2 = AF, AC BD. (36) 

We note that the weights of the two configurations shown in Fig. 7 are 
A 14 B 8 F 2 and A 15 B 6 F 3 respectively. If B 2 = AF, these make asymptotically 
equal contribution A l2 B 12 to X. 

Also, as A 2 F/B = AB, it follows that the weight increases by the same 
factor AB in going from configuration of type C\ to C 2 , and from C 3 to C 4 ( 
fig. 4). Similarly, the weight changes by the same factor BD/A = C, when 
going from C5 to Cq, as in going from Cj to Cg. Again, as tends to 
zero for large r, we have justified the asumption made earlier about such 
extentions being unfavorable. 

6 Critical exponents for b > 3 

If we start with the shortest polymer class of configurations that contribute to 
A^' +l ) having just weight A^, then the first extension is by a corner vertex 
at the boundary of the triangle, and the resulting weight is AB. For any 
subsequent extention, using either extension of type C\ — > C 2 , or C 3 — ■> C 4 , 
the weight increases by a factor AB. After n such extentions, the weight of 
the configuration is A n B n . For weight to be maximal, we want n to be as 
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Figure 8: Changing a dominant configuration for A^ r+1 ^ (middle figure) to 
one for _g( r+1 ) (left figure) or D^ r+1 ^ (right figure) by local change near a 
corner. The polymer is connected to another corner vertex of the (r + l)-th 
order triangle (not shown in figure). 



large as possible. It is easy to check that the largest value of n allowed ( to 
be denoted by j3 here) is 

= (b 2 - l)/4, if b is odd, (37) 
= b 2 /A } if b is even. (38) 

For example, both the configurations shown in Fig. 7, have f3 = 12. From 
the argument given above, any further extension of the polymer will be un- 
favorable. Hence such configurations correspond to the maximal term. If the 
number of such configurations is d(b), we have 

X = d{b)A p B p (39) 

It is difficult to determine d(b) explicitly for a general value of b. Of 
course, one can determine it for small b. This number would be expected to 
increase as exp(b 2 ) for large b. Fortunately, its precise numerical value does 
not matter for determining the exponents v and i/j. 

We can similarly determine y and Z. To modify the dominant configu- 
ration of A( r+1 ) to that for B^ r+1 \ one notes that as the former corresponds 
to a maximally extended polymer, one can find specific configurations such 
that the polymer reaches upto the r-th order triangle neighboring the corner 
which we want to reach. Then only changing the configuration in two such 
triangles is enough (fig. 8). This gives us 

y = B 2 /A (40) 
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Similarly, to get configuration of type /}( r+1 ) from X, we need only add an 
A-type vertex at that corner (fig. 8). This gives 

Z = A. (41) 

Putting all these together, we see that at the critical point, the recursion 
equations for and B^ are given by 

(\ogA^\_( (3 (3 \ / logAM \ / d \ (A0 , 
{ \ogB^ )~{f3-l (3 + 2 ) \ logflM ) + {c 2 ) {A2) 

For the 2x2 matrix given above, the equation determining the eigenvalues 

is 

A 2 -203+1) + 3/3 = 0, (43) 
whoes roots are 



A± = f3 + 1 ± ^f[3 2 - (3 + 1 (44) 

If 5 = x c — x is very small, then it increases by a factor A + under iteration, 
till it becomes sufficiently large so that linear analysis near the fixed point is 
no longer applicable. The number of iterations r max required for the deviation 
to become 0(1) is given by 

\ r ™*5 ~ 0(1) (45) 

whence we get 

r max = log(l/5)/logA + (46) 
the linear size of polymer varies as b rmax ~ (l/^)" with 

l/u = \og b \ + . (47) 

It is easy to check that for large b, v has an expansion of the form 

\lv — DJl — — — + higher order terms.) (48) 
olog(tr/2) 

For <Jc 1, logA^ increases as X r _. The maximum contribution to G(x) 
in equation (XX) comes from the term r = r max , hence we can replace the 
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sum by the largest term, and write logG(x) ~ logX^' max ^ ~ yy-ax^ ^his gives 
logG(x) ~ with 

(49) 

For large 6, this varies as log(3/2)/logb, and tends to zero as b tends to 
infinity. 

If G(x) varies as exp[(x c — x)~ a ] for x near x c , it is easy to see that 
coefficient of x n in the Taylor expansion of G(x) varies as x~ n exp(bn^) where 
b is some constant, and ip — a/(l + a). 

We have listed the numerical values of the critical exponents 1/V and ■?/>, 
along with the fractal dimension Db for some representative values of b in 
Table I. 



7 Discussion 

The critical exponent v does not tend to the two-dimensional value as b 
tends to infinity. This is not very surprizing, and a similar behavior has 
been encountered before in the case of the susceptibility exponent for linear 
polymers. Basically, there is a crossover from the two-dimensional euclidean 
value to the fractal value. For polymers with n monomers, with linear size 
<C b ( i.e. n <C b 1 ^, where v is the 2-dimensional value ), the space looks 
euclidean, and their mean size will be similar to that of polymers in regular 
2-dimensional space. However, polymers with n ^> b 1 ^ feel the constrictions 
of the corners strongly, and try to avoid them, and become more compact. 
Their average size is given by n v , with critical exponent v dependent on b. 
There is no exponent 9 that we can define for any b > 3, because of the 
presence of the essential singularity. 

One can define the chemical distance exponent z for our problem, just as 
we do for the euclidean problem. We take two sites on the branched polymer 
at a distance £ as measured along bonds of the polymer. If the average 
euclidean distance between these points is f(£), we define the exponent z 
by the relation f(£) ~ i 1 ^, for 1 £ <^ n 1 ^. We have not been able to 
calculate z for different values of b. 

We note that the logarithm of the number of configurations of polymer 
acts like the entropy. The non-translationally invariant fractal lattice pro- 
vides a deterministic model for the inhomogenous environment encountered 
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by the polymer in a random environment. In this case, it would seem more 
reasonable to average the logarithm of the number of configurations of rooted 
branched polymers over different positions of the root, as that would corre- 
spond to averaging the free energy of the polymer over different positions 
of the polymer in space. Such averages have been calculated only very re- 
cently for linear polymers on fractals [18]. It would be interesting to see 
if the stretched exponentional form for branched polymers is seen also in 
these "quenched" averages. Another open problem is the calculation of ex- 
ponents characterizing the collapse transition of self-attracting polymers on 
these fractals. It is hoped that future works will throw some light on these 
questions. 



Table I 

The fractal dimension D b , and critical exponents 1/V and ip for 

some values of b 



b 


D b 


l/i/ 


i, 


3 


1.63093 


1.41484 


0.13250 


4 


1.66096 


1.55263 


0.13381 


5 


1.68261 


1.57268 


0.12429 


7 


1.71241 


1.64448 


0.10702 


10 


1.74036 


1.70342 


0.09154 


15 


1.76787 


1.74406 


0.07825 


25 


1.79685 


1.78466 


0.06568 


50 


1.82788 


1.82292 


0.05375 


100 


1.85165 


1.84951 


0.04543 
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